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by 

Richard  L„  Liboff 

ABSTRACT 
In  the  classical  Chapman-Ens kog  theory  of  a  fully  ionized 
gas,  the  transport  coefficients  evolve  as  functions  of  certain 
weighted  total  scattering  cross-sections,,   If  a  Coulomb  potential 
is  assumed  between  scattering  particles,  these  scattering  integrals 
diverge,,   A  common  procedure  to  circumvent  this  difficulty  is  to 
integrate  the  Coulomb  results  up  to  the  Debye  cut-off „   A  more 
appropriate  procedure  involves  the  use  of  the  shielded  Coulomb 
potential,  in  which  case  the  scattering  integrals  converge,  and 
there  is  no  need  to  introduce  a  cut-offo   These  integrals  are 
calculated  using  the  shielded  Coulomb  potential  and  the  results 
are  compared  with  those  of  the  standard  "Coulomb  plus  Debye  cut-off" 
calculations.   To  the  "dominant  log  term"  the  two  procedures  give 
identical  results.   To  terms  of  the  next  order  of  accuracy  the 
difference  is  seen  to  remain  small „ 
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10      Introduction 

rn 

In  the  classical  method  of  Chapman  and  Enskog    for  the 
evaluation  of  the  transport  coefficients  of  a  fully  ionized  gas 
one  starts  with  the  Boltzmann  equation,,   Through  an  expansion  of 
the  distribution  function  together  with  solutions  to  the  subsequent 
linear  integral  equations  approximations  to  the  transport  coef- 
ficients are  generated,, 

These  approximations  to  the  transport  coefficients  entail  the 
evaluation  of  certain  weighted  total  scattering  cross-sections,.   If 
the  interaction  potential  between  the  ions  is  assumed  to  be  the 
Coulomb  potential,  then  these  weighted  cross-sections  diverge,, 
The  divergence  of  such  cross-section  integrals  is  well  known  in 
scattering  theory,  and  is  due  to  the  long  range  nature  of  the 
Coulomb  force o 

In  order  to  circumvent  these  difficulties,  a  common  procedure 
is  to  introduce  a  cut-off  parameter  in  the  limits  of  the  integrals 
instead  of  allowing  them  to  range  over  all  possible  distances 
between  the  interacting  particles,,   The  result  is  finite,,   However, 
the  choice  of  the  cut-off  value  is  rather  arbitrary  since  it  is 
based  solely  on  physical  arguments  which  are  dictated  by  the  nature 

of  the  particular  problem,, 

r  pi 

For  example,  Marshak, L  J  in  calculating  these  terras  for  star 

gases,  uses  a  quantum  mechanical  argument  to  obtain  a  cut-off 

value,  indicating  that  the  uncertainty  principle  would  prohibit 

the  use  of  a  classical  scattering  evaluation  for  scattering  angles 

beneath  some  minimum  angle,  which  he  chooses  as  the  cut-off  o   In 

Cut-off  in  angle  and  in  impact  parameter  shall  be  used  interchange- 
ably, due  to  the  identity,  sds  =  o~(cos  9)d  cos  9,  where  s  =  impact 
parameter,  o~  =  cross-section,  ©  =  angle  of  scatter,, 


such  a  calculation  one  fails  to  account  for  ions  scattered 

through  this  "diffraction"  region. 

r  71 
Landshoff    uses  a  similar  argument  in  determining  the 

cut-off;  however,  the  uncertainty  principle  in  this  instance  is 
used  to  discard  "head-on"  collisions,  so  that  the  scatterings 
through  angles,  9  #  %   are  excluded,  (opposite  to  Marshak's  use  of 
the  uncertainty  principle,  which  precludes  scatterings  through 
9^0).   To  remove  the  more  troublesome  scatterings,  9^0, 
Landshoff  exludes  collisions  which  last  longer  than  a  time  necessary 
for  the  ionized  gas  to  offer  a  "shielding  effect." 

Other  workers  in  the  field  have  suggested  that  the  inter-ionic 
distance  should  be  chosen  as  the  cut-off  while  still  another  choice 
is  the  "Debye  shielding  distance.   ^"   Roughly  speaking  this  dis- 
tance is  the  distance  within  an  ionized  gas  at  which  the  thermal 

2 

energy  is  balanced  by  the  electrostatic  energy. 

The  shielding  effect  was  first  mentioned  in  the  theory  of 
electrolytes  by  Debye  and  Huckel    .   They  showed  that  the  inter- 
action potential  between  the  ions  in  an  electrolyte  is  not  the 
Coulomb  potential  ~  —  ,  but  rather  a  shielded  Coulomb  potential 
of  the  form,  ~—  e~r'  ,  where  h  is  the  effective  shielding  distance, 
commonly  referred  to  as  the  "Debye  shielding  distance." 

More  recently  it  has  been  shown  that  a  similar  argument  holds 
for  an  ionized  gas.   If  one  assumes  such  a  modified  Coulomb  poten- 
tial rather  than  the  Coulomb  potential  in  the  above  method  for 
calculation  of  the  transport  coefficients,  the  related  divergences 


2  r^i 

It  is  interesting  to  note  that  Langmuir  ^  has  defined  a  plasma 
to  be  an  ionized  gas  for  which  the  Debye  length  is  small  compared 
to  other  macroscopic  lengths  of  interest. 

-  5  - 


in  the  weighted  total  cross-sections  disappear,,   For  those  tempera- 
tures and  densities  of  interest,  as  is  shown  later,  the  Debye 
distance  is  much  greater  than  the  inter-particle  distance,  which 
implies  that  a  binary  collision  approach  is  not  appropriate,,   How- 
ever, it  can  be  shown  that  calculations,  resulting  from  the  more 
appropriate  Pokker-Planck  equation  (derived  independently  of  the 
Boltzmann  equation)  agree  with  those  of  the  Boltzmann  equation  to 
the  "dominant  log  term„'^^^^ 

We  shall  concern  ourselves  in  this  paper  with  justifying  the 
usual  procedure  which  assumes  a  Coulomb  potential  together  with  a 
Debye  cut-off  in  the  evaluation  of  the  transport  coef f icients0   To 
this  end,  the  relevant  coefficients  evaluated  using  the  conven- 
tional cut-off,  will  be  compared  with  those  which  are  obtained  by 
assuming  a  shielded  Coulomb  potential,,   This  should  be  considered 
a  first  step  in  the  direction  of  obtaining  consistent  values  of 
the  transport  coefficients  since  the  assumption  of  a  spherically 

symmetric  shielding  cloud  is  only  an  approximation  to  a  more 

[91 
complicated  form  due  to  inertial  causes,,      Secondly  the  adopted 

potential  is  assumed  to  remain  unmodified  during  collisions,, 

The  author  wishes  to  acknowledge  his  grateful  appreciation  to 
Harold  Grad  and  Marian  H„  Rose  for  their  valuable  assistance  in 
preparing  this  report,. 


2 .   The  Transport  Coefficients 

The  following  formulas  for  the  various  transport  coefficients 
have  been  derived  in  the  usual  Chapman  Enskog  manner. 

For  a  two  fluid  gas  composed  of  ions  and  electrons,  so  that 
m  «  m. ,  it  can  be  shown  that  the  "second  approximation"  to  the 
coefficient  of  viscosity  (i,  reduces  to  that  of  a  single  gas  com- 
posed of  only  ions.   The  related  formula  appears  as 

(2.1)  ti  =  5kT/8  _fi_j2)(2). 

The  interaction  between  the  molecules  appears  in  the  calculation 
of  Xi-i   »  which  is  defined,  in  general,  by 

co         p 

(2.2)  SLU)lr)    -  /5E    f     e"V     V2T+2A(i)dV 

O 

(2.3)  A(^}     =      (   (l-cose0)sds 

O 

2       12 
where  V  is  a  non-dimensional  velocity,  given  by  V  kT  =  -p-mv  , 

9  is  the  angle  of  scatter,  (cf„  Appendix  1)  and  s  is  the  impact 

parameter,, 

The  viscosity  coefficient  jx  determines  the  thermal  conductivity 
coefficient  \   through  the  relation  \=2„7C^»   C   is  the  specific 
heat  at  constant  volume „ 

The  second  approximation  to  the  coefficient  of  diffusion  D. 
(subscripts  ie  indicating  diffusion  of  an  ion-electron  gas)  is 
written 


D   =    3kTM 
ie  "  l6nXlUJ(l) 


7  - 


where  n  is  the  number  density  of  the  particles  and  M  is  the 

reduced  mass  of  m,  and  m  „ 

i      e 

The  related  second  approximation  for  the  electrical  conduc- 
tivity  ©,  is  given  by 


ninen  ,  P 


2 

e  n  n 

— >    e 


kTn 


p  is  the  mass  densityc 


—  Die»   for  ^me/mi^  <K   1< 
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3o    The  Shielded  Coulomb  Potential 

The  expression  for  the  transport  coefficients  given  in  the 
preceding  section  show  that  the  problem  is  one  of  evaluating  the 
J\  integrals,,   In  the  introduction  we  have  already  stated  that 

these  integrals  diverge  if  a  Coulomb  potential  is  assumed  between 
particles.   However,  the  introduction  of  free  charge  (positive  or 
negative)  in  an  otherwise  neutral  plasma  will  induce  a  charge 
separation  which,  in  turn,  will  modify  the  potential  distribution 
about  the  charge „   This  modified  potential  takes  the  form  of  an 
attenuated  Coulomb  potential 9    the  attenuation  being  due  to  the 
shielding  effect  of  the  induced  charge  distribution,, 

The  equations  of  equilibrium  with  one  fixed  point  charge  in 
an  ionized  gas  are  written  (esu) 

(3d)     VPr  "  %E 

(3o2)    V°E  =  Mq  +  e16(r)) 

(3.3)    Pr  =  ^kTp, 
r 

E  is  the  electric  field  vectorn   The  subscript  r  denotes  the  r-th 

component  of  the  gas  (e„g0  for  an  ion-electron  gas  r  =  i  and  e). 

•6(r)  is  the  Dirac  delta  function,,   q  is  seen  to  be  the  "induced" 

charge  distribution  which  is  generated  by  the  fixed  point  charge  £-.<, 

Taking  the  gradient  of  (3°3)  gives  I 

kT 


(3.1U    VPr  "  rV^r 


r 
r 


(T   is  assumed  constant) 
r 


e  q 
r^r 


r 


Substituting   (3.1)    in   (3.10    gives 

(3.5)  7%  =    w^  E  . 

r 

Summation  over  r  yields: 

(3.6)  Zvqr  =  Vq    -  e  ^_     M 

of  which   the   divergence   is 

(3.7)  V2q  =  (V»E)Z^    +    E-Z^    ■ 

r  r 

We  now  assume  a  linear  theory  and  drop  the  second  term  in  the 
above  eauation.   The  remaining  equation,  together  with  (3.2) 
yields 

e  q 

2 

■ 

'r 


(3.8)   V2q  -  ^(l^Kq  +  ^sw) 


or 


(3.9)  v2q  -  \  (q  +  ex  6(r)), 

h 

(3.10)  h2  =  (Z_r_r)  sZ_^   o 

r  r 

The    solution  is    shown  to   be," 

en        -r/h 

(3.11)  q  =   -  — i-  ^— - 

Ipih2        r 

which  is  seen  to  have  the  correct  normalization,  viz0    qdf  =  -e. , 

oo 

If  one   introduces   a  potential   function  0   such  that. 


1  .9 


(3.12)      E  =  -70   , 


<2£  O  O      a" 

This  follows  from  the  well-known  identity  (V  °k  )     =  -Uac6(r) 

-  10  - 


then  the  related  Poisson's  equation  appears  as, 

(3.13)     V20  =  -Iptq  . 

If  q  Is  taken  to  be  the  total  charge  distribution,  (i„e<,, 
-r/h 


e,e 
[pik^r 

■r/h 


e.,6(r)  -  — - — 5 — ■  ),  then  the  solution  to  the  above  equation  is 
1         )iirV*T»  / 


e 


(3.110      0  =  e1 

and  the  corresponding  Coulomb  potential  is 

(3.15)      0     =  e-j/r  o 

As  already  mentioned  in  the  introduction,  the  length  h  was 
first  discovered  by  Debye  and  Huckel    in  investigations  concern- 
ing the  effective  fields  in  an  electrolyte „   The  concept  was  sub- 
sequently extended  to  plasmas  by  Langmuir,  Spitzer  and  others. 
The  form  (eq,  3.10)  derived  in  this  section  for  the  Debye  distance 

agrees  with  the  conventional  expression. 

1  /~T 
Equation  (3.10)  shows  that  h  ~  — /  —  „   This  is  in  accord 

with  the  interpretation  of  the  Debye  distance  as  being  that  length 

at  which  the  electrostatic  forces  are  balanced  by  the  kinetic 

forces,,   In  an  ionized  gas,  the  enormous  temperatures  present  are 

constantly  tending  to  separate  the  neutral  distribution,  while  the 

charges  tend  to  maintain  an  electrostatic  balance „   The  larger  the 

number  density  n,  the  less  the  mean  distance  between  particleSj, 

and  consequently  the  greater  the  electrostatic  forces  tending  to 

keep  electrical  neutrality,, 

The  ensuing  calculation  now  assumes  the  interaction  potential 


11  - 


between  two  scattering  particles  to  be  of  the  modified  form.   In 
essence  a  new  Boltzmann  equation  is  postulated  in  which  an  "effective" 
potential  is  substituted  for  the  usual  Coulomb  potential  between 
individual  particles. 

The  simple  rigid  sphere  model  of  a  gas  gives  for  the  mean 
free  path' ,  £, 

(3.16)  K   ~  -r- 

where,  for  an  ionized  gas,  the  molecular  diameter^  is  defined  by 
the  relations 

(3.17)  —  =  kT 

^  is  roughly  the  mean  distance  of  closest  approach. 
We  now  consider  the  product 

2 

(3.18)  AZ     =    e  a   =  ^~  =  h2 

kT^  n    ne 

which  shows  the  Debye  distance  to  be  the  geometric  mean  between 
the  particle  diameter^  and  the  mean  free  path  £„   If  we  introduce 
the  two  non-dimensional  quantities  £,  and  ^  ,  ass 


(3.19) 

«  -  A  A 

and 

(3.20) 

Z     =     K/h 

then  it  follows  from  (3.19)  that 
(3.21)     5~  =  l, 


For  a  plasma,  it  is  more  correct  to  take  the  mean  free  path  to  be 
the  above  value  divided  by  log  £(cf.  eq.  20),  also  [£]. 
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Ij..    The  A(  ^    integrals 

As  has  already  been  stated  In  part  2  of  this  report,  the 
calculation  of  the  transport  coefficients  involves  principally, 
knowledge  of  the  /v    integrals.   In  calculating  these  integrals, 
it  is  convenient  tc  introduce  a  length  L,  intermediate  between  the 
particle  diameter^  and  the  Debye  length  h,  i.e., 
(I4..I)      /\  «  L  «  h  . 

The  length  L  has  the  following  properties:  those  interactions 
which  occur  within  a  sphere  of  radius  L  are  governed  by  a  Coulomb 
potential.   For  interactions  with  impact  parameter  greater  than  L, 
only  grazing  collisions  are  important,  and  further  these  inter- 
actions are  governed  by  a  screened  Coulomb  field.   In  essence,  we 
are  assuming  a  new  form  for  the  potential,  namely, 

(lj.,2)  0    o*   i  r  <  L 

and 

(I1..3)  0  ~  expir/h!  r  >L 

which  is  seen  to  be  an  accurate  representation  of  the  modified 

Coulomb  potential  for  proper  values  of  L. 

A  ( 2.) 
Correspondingly  in  the  /v    integrals,  which  are  defined  by, 

00 

(k.k)  A{Ji)  =  J  sdsd  -  cos2e), 

o 
the  integration  is  divided  into  two  regions ,  so  that, 

L     oo 

(1+.5)    A(  }    =  J   +  J       • 

o     L 
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In  the  first  integral  the  exact  Coulomb  expression  is  used 
for  the  scattering  angle  9;  in  the  "far  region"  (s  >  L)  the  form 
for  0  corresponding  to  the  potential  (l|.3)  is  employed  together 
with  the  assumption  that  in  this  'far'  region  the  high  energy, 
small  angle  scatterings  are  dominant. 

As  is  well  known,  the  /\    integrals  diverge  logarithmically 
when  the  Coulomb  potential  is  assumed  throughout „   However  if  the 
modified  Coulomb  potential  is  used,  the  integrations  converge,, 

Before  evaluating  /\.        we  introduce  two  non-dimensional 
parameters  a,  and  o~,  through  the  equations, 

(kot>)  a   =  {£-)/?4-   , 

and 

( I4--7 )      «"  "  s/h, 

whence  we  see  that  small  a  is  indicative  of  high  energy  scatter- 
ing,  o~  is  merely  the  impact  parameter  measured  in  units  of  the 
Debye  length  h. 

In  the  expansion  of  the  first  integral  of  (k°5)    in  powers  of 
C,  the  leading  term  appears  as  (ah)   log(L/^)«  The  corresponding 

leading  term  in  the  expansion  of  the  second  integral  appears  as 

2 
(ah)   log  (h/L)»   Addition  of  the  two  series  immediately  indicates 

that  the  dominant  log  term  is  independent  of  the  choice  of  separa- 
tion parameter  L„  ((ah)2log(L/ ^  )  +  (ah)2log(h/D  =  (ah)2log(l/C)) 
The  next  order  of  accuracy  is  also  shown  to  be  independent  of  L0 

The  angle  of  scatter  9,  for  the  Coulomb  interaction  potential 
is  given  by  the  equation, 


-  lb 


(14..8)       csc2(|)  =  1  +  (2o/a)2. 

In  Appendix  I,  9  Is  calculated  for  the  modified  Coulomb 
potential  and  the  relation 

(l|.9)      ©  =  a  ^(o-  ) 

is  derived.  1L.    is  the  modified  Bessel  function  of  the  second  kind, 
first  order.   If  the  Coulomb  angle  of  scatter  (eq.  [4.08)  is  expanded 
about  a  =  0,  one  obtains  for  related  small  angles  of  scatter, 
©  ~  a/cr   ,   Likewise  if  K..  is  developed  in  eq.  (l\.t>9)    it  is  seen 
that  ©  ^  a/cr  for  cr~0,  so  that  the  result  for  the  modified 
Coulomb  potential  reduces  to  the  Coulomb  result  for  small  o~  .   In 
other  words,  far  within  the  shielding  cloud  (o~=  •?=  «  1),  the 
shielded  potential  reduces  to  the  pure  Coulomb  potential. 

We  are  now  prepared  to  calculate  the  /v    integrals.   As 
was  seen  in  part  II,  it  is  sufficient  to  consider  the  problem  for 
X  =  1,  and  JL-   2  only,  as  these  are  the  auantities  involved  in 
the  second  approximation  to  all  of  the  transport  coefficients, 
The  quantity  /\.    is  written 

CO  OD 

(J4..IO)     y\(2)  =  f  [sin2©]sds  =  h2  J  [sin2©]crdcr 

o  o 

or  in  accordance  with  (14-.5) 

<£  00 

(l4-.ll)     A(2)  =  h2  J  (sin2©c)ado-  +  h2  f  (sin2©c,  )0dcr 

^o   is  L  measured  in  units  of  Debye  length. 

The  subscripts  c  and  c>  refer  respectively  to  the  Coulomb 
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and  shielded  Coulomb  formulation.   Substitution  of  (^08)  and 

(2) 
c      cr 

•5; 


(ko9)  for  the  angles  9   and  9  .  gives  for  j\   '  , 


odo- 


which  upon  integration  yields. 


(U.13)        /^2)   =  i^|i-{iog(i+(2^/a)2)    -      (rff%2  +  2JK^K  +tfK   ): 

*    L  l+(2<£/a)^      -1  °    ° 


^v2} 


The  Bessel  functions  K  and  K,  are  evaluated  at  cr  =  <£ «   By 
definition,  o£  «  1,  and  furthermore  we  shall  restrict  the  calcula- 
tions  to  the  region  (2£/a)   »  1  (cf„  part  V  for  regions  of 
application) „   Using  these  simplifications  in  the  expansions  of 
the  Bessel  functions  and  log  term  in  (lj.«l3)   A.    becomes, 

(U-X1+)    A(2)  =  (ah)2  |ln(|)  -  1  =  Y  +  In  2J  . 

This  is  to  be  compared  with  the  expression  obtained  by  integrat- 
ing the  pure  Coulomb  formulation  up  to  the  Debye  distance  (^  =  lK 
Calling  this  form  J\\      ,   we  have, 

(14-.  15)   Aj.2)  =  (ah)2  |ln(|)  -  l]. 

These  last  two  equations  show  that  to  within  the  log  term 
the  Coulomb-with-cut-off  and  screened  Coulomb  formulations  give 
identical  results  for  the  Jy         integrals „ 

Recalling  the  defining  equations  for  a  (cf„  eq0  lj.o6)  and  the 
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variable  4  =  A /h  (cf«  ec*0  3.19),  it  is  easily  shown  that  a  may 
be  rewritten: 

(IJ..16)  a     =     4/2V2 

p  rr\XT 

where  V   is   the   non-dimensionalized  velocity  V     =  — g- /KPa 

This   expression  for   a   allows   eq's,    (^.ll;)    and   (l4-0l5)    to  be 
written  in  the   more   familiar   forms 


(M7)  AU)    =    (^h)2  (m(|)    +  2   ln(2Y)    -   .881*.} 

(1^.18)        A[c2)  =  <ah)2  {ln(£}  +  2  ln(2V)  "  ^   • 

1=1 

Prom  the  definition  (l+.lj.),  the  A(1)  integral  is 

OD 

(1|.19)    A(l)  =     (1-cos  0)sds 

o 

or,  more  conveniently, 

(iu2o)       A(1)=2J    sin2 1  sds° 

o 

Using  the  same  argument  as  for  A   ,  the  domain  of  integration 
of  A    is  divided  into  two  regions  by  the  length  L  giving 

L  „oo 

2 


(U..21)      fi1]   =2  f  sin2  §  sds  +  2    sin2  | 

o  L 


In  the  "far  region",  s  >  L,  the  approximation  is  again  employed 
that  only  scatterings  of  high  energy  and  small  angle  are  relevant, 
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2  9  9  2 

so  that  sin  ^     is  replaced  by  (^-)  „   Thus  the  second  integral  in 
c-  ,00  '   ^ 

1  f   2 

the  above  form  becomes  ■*  9  sds,  or  one-half  the  related  equiva- 

a(2)    L 
lent  part  of  the  yy    integral. 

For  the  first  integral,  we  have  from  the  Coulomb  formulation 


(U-o22)      sin2  %     =  — X    p 

2     1  +  (2cr/a)^ 


so  that 


S6 


(lt.23)     2  h2  f  sin2  §  o-do-=  £(ah)2ln(l  +  C2^-)2)  c 

w 
O 

This  gives  the  well-known  result  that  to  within  the  log  term  in 

the  expansion  of  the  /\    and  J)}        integrals  we  have 

c         c 

However,  if  all  terms  of  0(1)  are  included,  then  the  modified 
Coulomb  formulation  yields 

(^.25)     A(D  =  U|i!  |ln(|,  .  l.T  +  ln2}, 

whereas  the  ordinary  Coulomb  formulation  with  Debye  cut-off  (<£  =  1) 
yields 


(lv.26)  A<1>=     l#{ln|j 


In  terms  of  the  Debye  distance  h,  and  the  ionic  diameter  £±  , 
these  latter  two  results  appear  as5 
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(14-.27)  AU)  =  ^-    |m(|)  +  2  ln(2V)  -  \   -  Y  +  In  2J 

(14..28)  A^1)=  ^^   (ln(|)  +  2  m(2V)]0 

For  convenience  we  recall  the  parallel  results  obtained  for  J\ 
viz. , 

(1^.29)  A(2)  =  (ah)2|ln(^)  +  2  ln(2V)  -  1  -  Y  +  In  2J- 

(14..30)  AJ.2)  =  (ah)2/ln(|)  +  2  ln(2V)  -  lj 

(We  note  that  -y  +  In  2  =  .116). 
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$<>      Conclusions 

Ao   Regions  of  application 

The  major  restriction  imposed,  in  deriving  eq0  (I4-0I7)  of 
the  previous  section,  was  that  the  quantity  (2^5/a)   was  to  be 
large  compared  to  unity,  and  that  ~£    was  to  be  small  compared 
to  unity0 

Prom  the  relation  a  =  4/2V  f  2£/<x   may  be  rewritten  I4.XV  /<?<> 
The  quantity  E,   is  the  "particle  diameter"  measured  in  units  of  the 
Debye  distance 0   Expressed  in  terms  of  the  charge  e,  the  number 
density  n,  and  the  temperature  T,  g  is  written, 


(5.D         5-/5    5- 


;2 

kT 


,3/2 


■a       ■?/?  »A 

The   constant   ratio   e   /k  '    '  has    the   numerical   value     ^  10°°      (in  esu), 

2  2         2 

so    that    the    condition    (2^/a)      »   1,    or    (l+i&V  /£)      »   1,    may  be 

rewritten,  fTn  /o^l)      »  1. 
\n1/2lcTb  / 

_;£ 

The  value  of  X    is  chosen  by  realizing  that  e   /r  represents 

l/r  to  within  an  error  of  1  percent  if  X  is  10°*  „   The  non-dimen- 
sional velocity  V  is  taken  to  be  ~0(1)„   We  may  now  write  the 
above  condition  in  the  form, 

(5o2)       T3/n  »  10=8(°K  -  cm)3. 

From  this  relation,  the  present  formulation  is  seen  to  hold  for 
a  wide  variety  of  temperatures  and  densities  of  interest,. 

Bo   Results  and  remaining  problems 

We  have  seen  that  when  the  Coulomb  potential  is  assumed  to 
be  the  interaction  potential  between  particles  of  a  gas,  a  "cut-off 
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procedure"  is  necessitated  in  the  calculation  of  the  transport 
coeff icientso   The  cut-off  parameter  is  introduced  in  the  evalua- 
tion of  the  -/V    integrals,,   These  integrals  are  similar  to  the 

total  cross-section  integrals  of  scattering  theory,,   They  differ 

o 
by  the  weighting  factor  (1  -  cos  9),   Such  integrals  normally 

converge  for  inverse  power  potentials  of  the  form  r   ,  for  all 

values  of  n  greater  than  1„   In  the  Coulomb  case  n  =  1,  the  j\. 

integrals  (for  all  JL)   are  seen  to  diverge  logarithmically,, 

When  an  extraneous  charge  is  introduced  into  a  neutral  plasma 
a  shielding  cloud  is  generated  about  it  which  attenuates  the 
original  Coulomb  field  exponent i ally „   If  this  "modified  Coulomb" 
potential  is  used  in  the  l\         integrals ,  they  of  course  converge. 

In  the  conventional  calculation  of  the  transport  coefficients 
of  an  ionized  gas,  the  Coulomb  potential,  together  with  the  Debye 
cut-off  is  used0   We  have  compared  the  results  so  obtained  with 
the  transport  coefficients  evaluated  using  the  more  relevant 
shielded  Coulomb  potential,. 

The  results  show  that  to  the  "dominant  log  term"  there  is  no 
difference  between  the  two  calculations „   To  terms  of  higher 
accuracy,  for  the  diffusion  coefficient  and  electrical  conduc- 
tivity, the  related  error  is  of  the  order  of  two  percent,  and  for 
the  viscosity  and  thermal  conductivity  the  error  is  of  the  order 
of  one  tenth  of  one  percent „ 

The  present  results  have  been  derived  on  the  assumption  of 
a  spherically  symmetric  form  for  the  modified  potential.   As  has 
already  been  mentioned,  this  assumption  must  be  altered,  for  due 
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to  the  motion  of  both  scattering  particles,  there  will  be  an 

[91 
inertial  change  in  the  shape  of  the  screening  cloud.      In 

addition  the  modification  of  the  charge  distribution  during 

collisions  must  also  be  accounted  for» 

In  the  light  of  these  changes,  it  is  seen  that  what  has  been 

presented  so  far   is  a  first  approximation  to  a  more  general  and 

accurate  description  of  the  mechanisms  involved,, 
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APPENDIX  I.  Classical  Scattering  Formulation  for  a  Screened 

Coulomb  Potential 

Classically,  (when  the  orbit  is  well  defined)  the  angle  of 
scatter  for  an  arbitrary  potential  0   is  given  by 


(Al.l)     9  =  7i  +  2 


u 

sdu 


r°_ 

o   l/l-(jz!(u)/E)-(su 


7 


where  E  is  the  energy  of  interaction,  s  is  the  impact  parameter, 
u  is  the  inverse  radius  between  the  two  scattering  particles,  and 
u  corresponds  to  the  distance  of  closest  approach,  which  coincides 
with  the  root  of  the  form  in  the  above  radical,. 

Defining  the  three  non-dimensional  terms,  f(u),  a,  and  x, 
through  the  equations, 


(A1.2) 

0     =     J2fQ   f(u), 

(A1.3) 

a  =  0Q/E, 

and 

(Al.lj.) 

x  =  U//uo' 

eq.  (Al.l)  becomes 

1 
(A1.5)      9  =  %   +  2  f  dx 


o  /[(l-af(uox))/(suo)2]- 


x2 


The  equation  determining  u  is  written, 
(A1.6)     1  -  af(uQ)  -(suQ)2  =  0, 


of  which  a  first  iteration  about  a  =  0  gives 
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(A1.7)     1  -  a  f(l/s)  =  (suo)20 

Substituting  into  (Al<>£)  the  form  for  u  indicated  by 
(Al<>7)  and  expanding  the  resulting  equation  for  9  about  a  =  0 
up  to  linear  terms  in  a  gives  for  0 

r-  [f(u  x)-f(l/s)]dx 

(A1.8)      0  =  a  -/p . 

J      (I-X2)3/2 
o 

Within  the  order  of  accuracy  of  this  calculation  it  is  con- 
sistent to  replace  u   in  f(u  x),  in  the  last  formula,  by  ( l/s ) „ 

00 

We  may  now  write  for  9 

r1  [f(x/s)-f(i/s)]dx 

(1-x2)3/2      ■ 


(A1.9)       9   =   a 


-j 


o 
Another  form  of  (Alo9)  which  is  helpful  in  the  integration  to 
follow  is 

1  1 

U1.10)    9  =  -a  f  f(Vs)dx  +  a  {     [f(x/s)-x2f(l/s)]dx 

o   /^?      Jo      (1"x2)3/ 


(Al.ll)      =  -a(t/2)f(l/s)  +  a  f  [f  (xA  l^j  tl/s)  3dx 

J       (l-x2)J/2 

The  shielded  Coulomb  potential  is  of  the  form 

(A1.12)     0   =  A  eg^"r/h)  , 

"For  a  similar  discussion,  cf.  Kinetic  Theory  of  Gases,  E.  H, 
Kennard,  p„  119,  McGraw-Hill,  1938„ 


2h 


the  non-dimensional  part  f  being  the  above,  less  the  coefficient  A.' 
The  function  f(x/s)  appears  as 

(A1.13)       f  =  (x/cr)  exp(-cr/x), 

(crh  =  s) 

which  when  substituted  into  (Al„10)  gives  for  9, 

1 


it  \  /a\    /   \  .  fa< 


(kl.lk)  ©  =  -<§M§)  exp(-<r)  +  (£) 


(x  exp(-o/x)-x  exp(-p-)) 


dx, 


Employing  the  substitution  x  =  cos  0,  followed  by  a  parts 
integration,  yields  for  (Al.ll|.) 

Tt/2 
(M.1S)    e=-aj  M.*  •*.<— «*• 

o 

The  desired  result  follows  immediately  from  the  integral 
formula 

it/2 

(A1.16)        secn0  exp(-crsec  0)&0  =   (-l)n'1K^n"1)  (cr  ), 
o 

for  the  (n-l)st  derivative  of  the  zero-order  modified  Bessel 

function  of  the  second  kind  K  .   (Exponents  without  parentheses 

are  algebraic.)   The  above  formula  together  with  the  relation 

K  (o-  )  =  -YL  (o-  )  reduces  (A1.15)  to 
o         1 


(A1.17)  9   =   aK^cr  ). 


'cfo  eq.  (3.1l4-)o 

,"00e"Xtdt 

This  formula  follows  from  the  relation  Kn 


^■Jj 
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The   expansion  for   K_  ( cr)    is 

oo      /_ /„x2i+l 


(A1.18)      K-^cr)   =   (|)   +Z     u(iil)l    {log(o-/2)-(l/2)[^(ifl)^(i+2)]]p 
where   the   \|/   functions    are   given  by 

(A1.19)  \|f(i+D    =   (1  +  |  +  ^  +    000   +  i)    -   y,    *(D    =   -Y. 
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APPENDIX  II,  Expansions  Involved  in  the  TV    integral 
The  7v    integral  was  written 


oo 


(A2.1)    A(2)  =  f  (1  -  cos26) 


sds 


oo 

J  ■ 


sin  9  sds  . 


The  integral  was  approximated  by 

L 
(A2.2)    _A(2)  = j 


J" 


2         I    2 
sin  9  sds  +    ©  ,  sds. 
c     J       c' 


.  2, 


In  the  first  integral,  the  Coulomb  formulation  for  sin  9   gives 


(A2o3) 


/  sln 


9  sds  -  — ?5 — 
c        2 


(2^)2 

iog(i+(— ) ;-  i+(2^)2 


and,  to  within  the  relevant  approximation  (2^/a)   »  1,  this  part 
of  the  integration  is  approximated  by 
L 


(A2.1U 


/ 


.  2n   ,     (ah) 
sm  9  sds  =    0 
c        2 


2  r 


2  1 


og(f ,  -  l] 


The  remaining  integration  over  9  t  is  written 


(A2.5)       ©c'  sds  =  (ah)2    I^(<r)o"d<r  . 
L  <£ 

This  integral  is  expanded  as 

oo 
(A2.6)    J  o-K^(o-)do-  =|X2  J42(£)  -  (1+  ^  )xfOC)J 

Using  the    identity 

(A2.7)  £2K{2      =      (^   +    2KQ)2, 


-  27 


allows   the    above    integral   to   be   rewritten 

oo 
(A208)  o-K^dC    =  \  |-£2K^     +   2£k]K0   +    X2Y?\  . 

<£> 

Expansion  of   the    three   terms   in   (A2„8)   yields, 

(A2.9)  X2K^  =    (1  +  -4"  ln  T"      +      •••    )2     =       l+0(X2ln^), 


(A2.10)  ^K2   =    (^2Y  +    X2    In  ^  +    ...)2  =0(Xk  In  ^-) , 

and   for    SoK-K , 

(A2.ll)         X\XQ  =   -l&r  +    ^ln  ^  +    ooo)(X     +  it  ln(-T}    +    °oo) 

=   -    (Y  +  In  ■§-)   +    0(^2   ln  ^0. 
These   expansions  permit   the   integral   (A2„8)    to  be   rewritten  as 


oo 


(A2.12)  trader     =      -|-Y-ln^-+    0(£2   In  <=§    ) 

Combining  the    two    integrals   in   (A2»2)    gives   for  Jv      , 

(A2.13)        Ai2)   =    (ah)2  [lnl2^)    -  |  -  |  -  Y  -   ln(4)+0(X2   In  f )j 

Retaining  only  terms  within  the  brackets  of  the  order  of  unity, 
yields  finally  for  A , 

=  (ah)2  |ln(|)  -  1  -  Y  +  In  2 

=  (ah)2[ln(|)  -  .881+1  . 
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